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^ ■ Abstract 



Using the unique nature of the average output state of an optimal signalling ensemble, 
we prové that for a special class of qudit unital channels, the HSW channel capacity 
is C = log 2 (oí) — min p S {£{p)), where d is the dimension of the qudit. The result is 
extended to products of the same class of unital qudit channels. Thus, the connection 
between the minimum von Neumann entropy at the channel output and the transmission 
rate for classical information over quantum channels extends beyond the qubit domain. 



O 

Oh. 1 Introduction 

H— > 

The Holevo-Schumacher-Westmoreland theorem telis us the asymptotic rate at which clas- 
sical information can be transmitted over a quantum channel £ per channel use is given by 
the maximum output Holevo quantity x across all possible signalling ensembles. 

•i-H . 

^ ■ C = max[ p . íPi} x({pí,Pí = £{Pi)}) 

Here x is the Holevo quantity of an ensemble {pi,p%}, defined as 



X 



5 (^L,PiPi) ~ Pi S (Pi 



where S is the von Neumann entropy. [] We call C the Holevo-Schumacher-Westmoreland 
(HSW) channel capacity. 

We call any input ensemble {pi , p{\ that achieves C an optimal ensemble. There may 
be several different optimal input ensembles which achieve the optimum HSW channel 
capacity C. However, it was shown in [Q] that the average channel output state of an 
optimal ensemble is a unique state for all optimal ensembles for that channel. That is, 



We shall use p to denote a density operator at the channel input, and p as the corresponding channel 
output density operator. 
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given a set of optimal input ensembles ip- , p\ j , <.p\ , p \ > , ■ ■ ■ , yp\ , p \ >, all 
achieving C, we define $W = £ P{^ P^^j ■ Then it has been shown we must have 

$(1) = $(2) = ... = $(A0_ 

The main idea of this paper is the unique nature of the output ensemble average state of an 
optimal signalling ensemble for a quantum channel £ telis us alot about C for that channel. 



2 Background Material 

2.1 Invariance of S and \ under unitary operators 

Consider any ensemble {pí,Pí}- Acting on each pi with the same unitary operator U yields 
a set of vàlid quantum states UpiU^ and the ensemble UpiU^j. Furthermore, each pi 
has the same eigenvalues as the corresponding UpiU^ . Since von Neumann entropy depends 
only on a density operators eigenvalues, we conclude S(pi) = s(lJpiU^. Furthermore, 

this implies the Holevo quantity x of the ensembles {pí,Pí} and jpj, UpiU^ is equal, since 



x({pí,u Pi uï}) = s(^Ç Pl u Pl u^ - Y,PiS(u Pi tf) 
= s(u(^Ç PiP \ tf\ - Ç piSiypirt) 

= S {^2píPí^ - ^PiSiPi) = x({pi,Pi})- 



(I) 



3 HSW Channel Capacity for single qubit unital channels 



As an example of the approach we shall be taking, we derive the HSW channel capacity 
for single qubit unital channels. This result was previously derived in [Q] by a different 
technique. 

We describe a single qubit density operator using the Bloch sphere representation. 

p = ~ (X + W p .ff) 
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The symbol a is the vector of 2 by 2 Pauli matrices 



u 



where a 7 
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The Bloch vector W is a real, three dimensional vector which has magnitude equal to one 
when representing a pure state density matrix, and magnitude less than one for a mixed 
(non-pure) density matrix. 

It was shown in |^] that the action of a single qubit unital channel £ on an input state 

'Wx 

p could be represented as p = £(p), where p has Bloch vector 
A x w x 



w 



w z 



vector 



X y w y 

\ z w z 



and p has Bloch 

Here the Àjt € [—1,1]. Using the unique nature of the average output 



state of an optimal signalling ensemble, we shall show the HS W channel capacity C is 
C = 1 — maxp S(£(p)). 



3.1 Achievability of Output Ensembles 

We say an ensemble {qj , (frj} at the channel output is achievable if there exists an input 
ensemble {qj , cpj} such that the {fj} are all vàlid density operators and £(fj) = 4>j Vj. 
Let us recali some properties of the Pauli matrices {cr^}. The obey the relations Gi(Jj = 
— OjOi for i ^ j and a^aj = I2 for i = j. Thus, we find aiajai = — o~j for i ^ j and 
OiO-jOi = ai for i = j. The are Hermitian, so a\ = I2 implies the are unitary, 
yielding a\ = a k . 

Let cebe an optimal input ensemble with corresponding output ensemble {pi , £{pi) = Pi} ■ 
Apply a Pauli operator u\- to all the density matrices in {pi , £ (pi) = p,j} , yielding 
an ensemble {pi , o^piO^}. We know the density operators {o~kpiO~\} are vàlid density 
operators because ak is a unitary operator, and hence acting with ak implements a change 
of basis at the channel output. The question we are interested in is whether the output 
ensemble {pi , ak Piaj,} is achievable. To answer this, we know for each pi, there is a vàlid 
input pi such that £(pi) = p%- Consider the following. 

akPia\ = a k £(pi) a\ = a k £ Q(^2 + u Xi o x + uJ Vi a y + uj Zl a z ^j a\ 
= a k {-zíh + Xxio Xi a x + \ y w yi a y + \ z uj Zi a- 



"We write Id for the d by d identity matrix. 
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= -{j2 + \xUxi<TkVx<r\ + ^y UJ y l (T kO'yCrl + Xz^z^k^zO'l 
Define S^i = if k = l, and 1 if k 7^ l. Note that cr^aia^ = (— l) 5fe,i 07. If has the 



Bloch vector 



(-1)**.- Cü x 



then the channel output of ifi is 



e{(pi) = \ {h + {—l) &k ' x X x uj x a x + (-l) 5 ^y X y u y a y + (-l) 5 ^ X z uj z a z 

= 2 + ^x L °x' J kO'xO'Í + \UyO k O y o\ + XzU)z<?k a z&k) = Vk^iVi) a \- = ^kPi a \- 

If we can show the (fi are vàlid density operators, then we have shown that the output 
ensemble {pi , a k pi a\} is achievable. In order for ifi to be a vàlid density operator, we must 
have the corresponding Bloch vector composed of three real entries, and the magnitude of 
the Bloch vector less than or equal to one. Since the pi are vàlid density operators, the 
three ui k are real, and obey uj x + u y + u> z < 1. Now (— l) 5fc '' for k, l = {x, y, z} is real and 



Wi ) + 



\Sk,x 



equal in magnitude to one. The magnitude of the Bloch vector for ipi is í (—1 

( {-lf k >y tüyY + ( (-lf*'* = u? x + u? y + u? z < 1, where the last inequality 

follows from our knowledge that the pi are vàlid density operators. Thus the ifi are vàlid 
density operators. We conclude that if there exists an optimal input ensemble ce, with 
corresponding output ensemble {pi , £{pi) = Pi] , then the ensemble jpj , o~ k pio\}^ is 
achievable, with corresponding input ensemble {pi , (fi}. Furthermore, the input ensemble 
{Pi , fi} is optimal, since a k is a unitary operator, and we showed in equation (|) that 
a unitary operator acting on an ensemble does not change the Holevo quantity of that 
ensemble. Since {pi , £{pi) = pi} attained the maximal Holevo quantity C at the channel 
output, the output ensemble ipi , a k pia]}^ also has a Holevo value of C. Thus {pi , ifi} is 
an optimal input ensemble. 

To summarize, we first chose a basis of operators Ei, in this case the identity I2 and the 
three Pauli operators {a x , a y , cr 2 }, in which to expand the density matrix p = J2i a iEi- 
Next, we found a set of unitary operators U/., in this case again the Pauli operators a k , such 
that the U k act on the Ei resulting in a multiplicative phase factor : U^Eil/j, = Kt k A Eí, 
where Kn.^ is a complex quantity. The unital nature of the qubit channel £ telis us that 
£{Ei) = XíEí Vi in the operator basis {Ei}. This leads to the commutation of the 
channel £ with the set of unitàries {U k } = { ± h , ±& x , ±a y , ±c 2 }. 

U k £ {Ei) U\ = UkhEiUi = \U k EiUl = \ « (M) E { = k (m) £ (E t ) 

( By linearity of quantum channels ) = £ (k^,í) E^j = £ (u k Ei . 
Since we have an expansion of p in terms of the Ei, using the linearity of quantum channels, 
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we conclude that 



U k £(p)UÍ = U k £ (~£ a ^J 4 = U k = 

(II) 

= I^oíS^EíUI) = S^^aiUkEiU^j = £ ( U k ^ Ç «<^J ^í) =s(ü k pUl) 

A Uk acting at the input is a basis change and hence UkpUl is a vàlid input density 
operator. Equation allows us to conclude that any Uk acting on the output states pi of 
an optimal ensemble oeyields an output ensemble {pi , UkfiiUl} which is achievable. The 
achievability of channel output ensembles generated by U k acting on the output ensemble 
of an optimal input ensemble will be a critical tool in extending our unital qubit channel 
analysis to the determination of the Holevo-Schumacher-Westmoreland channel capacity C 
for a special class of qudít unital channels. 



3.2 Using Symmetry Properties of Optimal Ensembles 



Consider a unital qubit channel with an optimal input ensemble {pí,Pí}, f\ average input 
state <ï> = J2i PiPi an d average output state $ = £ ($>). Let $ have Bloch vector V = 

V x ' 

Vy 

v z . 



and <1? have Bloch vector V 



~V X ~ 




' X x v x ' 




Vy 




Xy Vy 


. Choose one of the three {cx^} 


.Vz. 




,X z v z _ 





and apply this Ok to the output states pi to obtain a new output ensemble {pi , OkPio\} = 
{pi , p'j}. We know from our work above that the output ensemble {pi , <JkPi<j\} is achievable 
and optimal. The action of Ok on the output ensemble {pi , £(pi) = Pi} generates a 
corresponding transformation of the average output state of the optimal ensemble 



E 



Pi O k Pi o\ = Ok 



E 



Pi Pi o\ 



Ok^O 



By the invariance property shown in pj, we have = <è. Now $ has Bloch vector 



V 



~v x ' 




Ax v x 


Vy 




Xy Vy 


.Vz. 




_X z v z _ 



{x, y,z} , <ï> = implies 



1 \&k,x 



and 







"(" 




V x 


Bloch vector V = 




(" 




Vy 






u- 




V z . 


— l) 5fe -!' v y and 


V 


z - 


- ("1 





For k 



(III) 



The only way all three relationships in equation ( p| ) can be true V k = {x, y, z] is if 
0. The fact <3? has Bloch vector V 



leads to the conclusion that 



3 That such an ensemble exists was shown in pj. 
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<£» = g ( 1% + V • o ) = i/2 for all optimal ensembles. 



A second way to see that$ = \h is via Schur's lemmaQ. Consider the group TL composed 
of the eight operations { ±Í2 , ±o~ x , ±a y , ±a z }. A necessary and sufficient condition for a 

finite group Q to be irreducible is if the relation X^eç Trace[g] = 1 is true[Q]. Here 
||Ç|| is the order of the group Q. Noting our group TL above is finite, and computing the 
sum with our group TL, we find TL is irreducible. 

Schur's lemma states that if a group Q is irreducible and has a d dimensional representation 
r {G) in which each representation element T{g) commutes with a d by d matrix M \/g £ Ç, 
then M is proportional to IdP]. The fact that we found a^a^ = $ VA: S {x,y,z}, 
together with the same trivial result for I2, implies that all elements of TL commute with $ 
and thus $ oc I2 ■ The trace condition Trace (^j = 1 leads us to conclude <£> = \ I2 ■ 

Having determined we can now rewrite the Holevo-Schumacher-Westmoreland channel 
capacity C as C = log 2 (2) — J2i Pi $ (£ {Pi}) ■ To further simplify this result, we use two 
results from ||]. In their paper, Schumacher and Westmoreland worked with the relative 
entropy function, X?[jo||çí>] defined as Trace [p\og 2 {p) — / olog 2 ((/')]· Using V, they proved 
the following two results. 



I) The equal distance property of optimal ensembles. 
For any optimal ensemble ce, we have 



v 8{ Pi ) e (*) 



C Vi. 



(IV) 



II) The sufficiency of the maximal distance property. 

For any optimal ensemble oewith average input state $ = Y^íVíPí-, we have 
T> £(4>) £(&) < C for any input density matrix (j) 



(V) 



In both I) and II), $ = J^íPíPí an d C is the Holevo-Schumacher-Westmoreland channel 
capacity. For the case of qubit unital channels, we have found that every optimal ensemble 



cemust 



obey £ ( £, 



Pi Pi 



Looking at the relative entropy formula, we see that 

V £{4>) \l d = log 2 (d) - S{£{4>)), where S is the von Neumann entropy and (f) is any 
input density matrix. (For more details on this derivation, please see the appendices in Q.) 
Using the fact that for qubit unital channels we have found, for all optimal ensembles ce, 
that £(^J2íPíPí^) = è-^ 2 ' above two Schumacher and Westmoreland results become, in 
the qubit unital channel case, 
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r) 



1 - S{8{pi)) =C Vi implying S (£ (pi)) = S (£ (pj)) 



(VI) 



II') 



1 — 5 (£(</))) < C V input density matrices ep. 



(VII) 



We know that II') is achieved with equality when (p is any of the pi in the optimal ensemble oe. 
Thus I' ) and II') taken together yield 1-5 (£ (</»)) < 1 - S {£(pi)) or S (£ (ep)) > S(£(pi)), 
which, since (p can be any input density matrix, implies S (£(pi)) = min^ S(£((p)). 
Plugging this result into I') yields our final result for the Holevo-Schumacher-Westmoreland 
channel capacity for qubit unital channels. 



For qubit unital channels, the minimum channel output von Neumann entropy determines 
the Holevo-Schumacher-Westmoreland channel capacity C. 

3.3 Ensemble Achievability 

The achievability of a transformed output ensemble is a concept worth emphasizing. In 
our discussion of unital qubit channels, the reason why we could conclude the average 
output state of an optimal ensemble commuted with all eight members of our group H = 
{±Í2 , ±t x , ±a y , ±a z } was because, given an optimal ensemble ce, each of the eight output 
ensembles {pi , hpih^ 1 }, where h G H, was achievable. The existence of an optimal input 
ensemble {pi , (pi} which maps via the quantum channel £ to {pi , hpih^ 1 } is what allowed 
us to conclude the relationship hQh -1 = <3? was vàlid, and apply Schur's lemma. 

For a genèric group Ai acting on the channel output of an optimal ensemble oe, there will 
typically be mo G M such that jp^ , mopirriQ 1 ^ are not achievable ensembles. In these 

cases, we cannot conclude mQ^rriQ 1 = í> holds, where $ is the average output state of 
an optimal ensemble. Yet it was the fact that mo^rriQ 1 = <ï> holds Vm G M. that led us 
to apply Schur's Lemma and conclude $ oc Í2- The lack of achievability for one or more 
of the transformed output ensembles jpj, m pi m" 1 j prevents us from appealing to Schur's 
Lemma. An example of the limitations to determining HSW channel capacity which results 
from output ensemble non-achievability arises in the case of non-unital qubit channels. 

3.4 A Non-Unital Qubit Channel Example 

Our technique fails for non-unital qubit channels. The reason why is the lack of achievability 
of output ensembles generated by members of the Pauli group acting on an output optimal 



C = 1 — min^S (£(4>)) . 



7 



ensemble. For example, consider the non-unital linear qubit channel specified in the Ruskai- 
King-Swarez-Werner notation as {t x = t y = 0,t z = 0.2, X x = X y = 0, X z = 0.4}. This 

channel maps an input Bloch vector W to an output Bloch vector W as : 



W 



Wx 

w 



w. 











0.2 + 0Aw z 



W. 



By inspection, an optimal input ensemble is cewith pi^ = \ {h ± c 2 ), and corresponding 
output density matrices p\ = | (I 2 — 0.2<r z ) and = \ (h + 0.6a z ). Numerical analy- 
sis for this channel indicates the optimum output average state is <£> » \ ( I 2 + 0.2125 a z ). 
Since $ 7^ ^ I2 , we anticipate we will not be able to meet the conditions for the application 
of Schur's lemma. 

Consider applying the unitary operator a z to the output optimal ensemble {pi , £{pi) = 
Pi} determined in the previous paragraph. We obtain 



&z Pi Oz = dz Q (^2 - 0.2a z ^) ja\ = p\ and a z p 2 cr z = cr z Q (l 2 + 0.6<7 Z ) ^ a\ = 

Thus the output ensemble jpj , o z E {pi) o\ = cj 2 pjcj|| is identical to the output ensem- 
ble jpj , 6 (pi) = pi j, both being generated by the input ensemble ce. Thus the output 
ensemble jp» , a z £ (pi) a\ = a z pi <rj j is an achievable output ensemble. 

The application of a x or a y to {pi , £(pi) = Pi} however does not yield an achievable 
ensemble. To see why, consider applying a x to P2 = \ (h + 0.6cr z ), which since cr x a z a\. = 
—a z , yields the output density operator p~2 = \ (h — 0.6a z ). The corresponding input 

r o" 

density operator would have Bloch vector W = , which is not a vàlid qubit density 

L-2. 

operator, since \\W'\\ > 1. Since the output state a x p~2<j x can never be mapped to by a 
vàlid input qubit density operator, we cannot assume the relation a x Í a x = $ holds. 
Thus, we do not have the necessary Schur commutation requirement that g& = &g for all 
members g of the Pauli group {i/2, ±<7 X , ±o"y, ±<7z}, and hence cannot conclude = \l 2 , 
as we anticipated. 

As we shall develop in more detail below, working with qudits, if we can find a group Q 
such that we are assured all elements g G Q are unitary and acting on the output states of 
an optimal ensemble ceyield achievable ensembles V g G G, than we will be able to conclude 
the average output state of any optimal ensemble is $ = -j Id- From this conclusion, 
we can use the Schumacher-Westmoreland relative entropy results from equations |V],[V],|VÏ|, 



and VII to conclude the states in any input optimal ensemble must be a subset of those 
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input states which yield the minimum output von Neumann entropy. This in turn leads us 
to a HSW channel capacity C of 

C = log 2 (d) - min^S (S(4>)) 

for those qudit channels to which we can successfully apply Schur's lemma. We now proceed 
to determine the subset of qudit channels which meet the Schur's lemma requirements. 



4 Qudit Channels 

The HSW channel capacity result for unital qubit channels was previously proven in || by 
a method which did not generalize to the general qudit case (ie: for qudit dimension d > 2). 
The technique discussed in this paper does generalize to a special subclass of unital qudit 
channels. Before describing that generalization, we present some background material on 
qudits and qudit channels. 

4.1 Qudits 

A qudit is a system with d orthogonal pure states j = 0, 1, 2, • • • , d — 1. The generaliza- 
tion of the qubit Pauli operators a x and a z are the two operators X and Z, whose action on 
the states |j) is X\j) = \j + 1 (modd)} and Z\j) = ÍP \j). Here Cl = e~ . The extension 
of the qubit Bloch representation for a density matrix p to qudits is shown in appendix A 
to be 

P = \ E X a 

a,6e{0,l,2,-,d-l} 

The a a ^ are complex quantities. Define E a ^ = X a Z b . Note that Eq % q = Id- I n appendix 
A it is shown Trace(E a ^) = dò a fl 3b,o, where 5 is the Kronecker delta function. The trace 
condition Trace(p) = 1 al·lows us to conclude ao,o = 1- Let T denote the set of d 2 — 1 
elements a, b £ {0, 1, 2, • • • , d — 1} with the exception that a and b cannot both be zero. 
Then we can write the qudit density matrix p as p = \ (jd + J2(a,b)eT a a,b -^a,è) • A 
qudit quantum channel £ is a linear map. One can write such a map as a d 2 by d 2 complex 
matrix A4 taking the d 2 vector of coemcients a a b of p to the d 2 set of coefhcients à a t of 

p = £(p)- S 

If the qudit quantum channel £ is unital, meaning £{Id) = Id, then the first row and 
column of Ai are a one followed by d 2 — 1 zeros. Hence we can represent a qudit unital 
channel by a matrix N of d 2 — 1 by d 2 — 1 complex entries mapping the vector of d 2 — 1 

4 Our qudit matrix development in which we write £ as a d 2 by d 2 matrix closely follows work done in || 
for the unital qubit channel case. 
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coeficients au b \, with (a,b) G T, representing /) to the vector of d 2 — 1 coeficients 5( ai 5), 
with (a, 6) G T, representing p = £(p). The specific class of qudit channels we shall be 
interested in are those completely positive unital quantum channels for which J\í is diagonal. 
This class of channels is nonempty. For example, consider the channel corresponding to all 
zeros on the diagonal. This point channel maps all input density matrices to a single output 
density matrix p = \ld- Another member of the set of diagonal unital channels is the 
identity map, which maps any input density matrix to itself. This channel has all ones on 
the diagonal of the matrix Ai. 



The approach we take to determine the HSW channel capacity for this special class of 
diagonal unital qudit channels closely follows our unital qubit channel derivation above. 
Note the operators E a b are unitary. Using the commutation relation shown in appendix A, 
ZX = ÍIZX, where il = e d , we have 



E g ,h E a , b E h 



x 9 z h x a z h z- h x~ 9 = Çl ah X 9 X a Z h Z b Z~ h X~ 9 



(VIII) 



= n ah x 9 x a z b x~ 9 = ü ah n~ bg x 9 x a x~ 9 z b = n ah n~ b9 x a x 9 x~ 9 z b 
= n ah n- b9 x a z b = n ah ~ b9 E a)b . 

Define F a ^ c = £l c E ajb , where a, b, c G {0, 1, 2, • • • , d — 1}. Since O c and the E a ^ are unitary 
operators, F a;b)C is a unitary operator. The action of the F a ^ fi on a diagonal unital qudit 
channel output density operator p is 



(IX) 



E n 



h + ^1,r Oíq,r E q ,r J E a h 

(9,r)6T 



d [ Id + E x i> r a i' r Ea > b E w E l ;b 

(q,r)eT 



2( Id+ E ^«^1=^1^+ E a q , r n b9 ~ ar E q , r 
(q,r)eT ) \ \ (g,r)eT 



= £ E a>b - \l d + ^2 a q>r E q 

V V (5.r)eï 

Since the 6 c are unitary operators, we conclude that given any optimal input ensemble ce, 
the output ensemble @ a ,b,c obtained by applying F a<bíC to {pi , £{pï) = pï\ is achievable and 
@ a ,b,c has the optimal input ensemble jpi , = F a ^ c pi F\ b \ . Each of the (pi is a vàlid 
input density operator due to the fact that F ajbtC is a unitary operator, and is implementing 
a change of basis on pi . 

The set of operators {-F a ,fe,c} forms a group of order d 3 which we shall call Q. Recali 
our theorem for proving a finite group is reducible. |0| The group Q is reducible since 



Trace 



Faj 



equals zero when either a and b are non-zero, and 



Trace 



b.c 



equals 
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2 

= -jjdd = 1. Since Q is a reducible 



d when a = 6 = 0. Thus Trgjr EggQ Trace[g] 

group, we can apply Schur's lemma. For any optimal input ensemble ce, the channel out 
put ensemble jpj , F a ^ c £{p,i) F^ h \ is achievable and the corresponding input ensemble 

\pi , Fa,b,cPi -^Ibc} is optimal. From the uniqueness of the average output state $ for any 
optimal ensemble, we conclude that V a,b, c : F a> b íC <3? F^ b c = <ï> or F a ^, c $ = ^^a.fe.c- 
By Schur's lemma we obtain $ cx 7^. The trace condition telis us Trace (^J = 1, so we 
conclude <£> = \ld- 

This leads us to conclude that for the optimal input ensemble ce, the HSW channel capacity 
is C = log 2 (d) — J2i PiS{£{Pi))- Using the two relative entropy results from equations 
^V|JV|JVÏ|, and |yTJ, as we did in the qubit case, we obtain S{£(pi)) = min^ S (£(4>)) yielding 
the HSW channel capacity for diagonal unital qudit channels. 

C = log 2 (d) — min^S {£{(j))) . 



5 Products of Diagonal Unital Qudit Channels 

Consider the product of N diagonal unital qudit channels £^\ k = 1, • • • ,N. The tensor 
product channel is £ ® = £^ <S> £^ <£> ■ ■ ■ <8> £^ N \ Let the input qudit density operator p^> 
corresponding to the diagonal unital channel £^ be of dimension df~. Then d = Yl^=i dk 
is the dimension of the input qudit p® for the product channel £® . The basis elements for 
p® which we shall use are the tensor products of the individual E^f) . 

{>%} = 

where the ük and bk £ {0, 1, 2, . . . , dk — 1} and a and b E {0, 1, 2, . . . , d — 1}. 

The basis elements E® b are not necessarily constructed using the d dimensional qudit op- 
erators X and Z described in appendix A. As a result, we must prové several properties 
for the set {-E^f;,} before we proceed with the HSW channel capacity analysis for product 
channels. 



5.1 The mapping {oi, a 2 , (13, ■ ■ ■ , ün-i, o-n} {«}• 

The set of possible coefficients {ai, a 2 , 03, • • • , «at-i, cln} and the set of possible {a®} both 
have d elements, where d = n|=i dk- Here the {a^} £ {0, 1, 2, • • • , dk — 1} and {a} G 
{0, 1, 2, • • • , d}. There are many bijective mappings between these two sets, and it is useful 
to have one particular map in mind as we proceed. The one we shall use is presented in the 
table below. 
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a N 






































ajv-i 


































































«3 






































0,2 
















1 


1 


1 




1 


2 


2 




a\ 





1 


2 




di - 1 





1 


2 




di - 1 





1 




a 





1 


2 




di - 1 


di 


di+1 


di +2 




2di - 1 


2di 


2di + 1 





and so on. Below, we associate an E a b with the tensor product ^-E' ai bl C5> 6 2 ® ' ' ' ® E a^ 6 
by using this mapping twice, once for the association {01, 02, 03, • • • , ajv-i, o>n} < í=^ {a®} 
amd again for {61, Ò2, 63, • • • , &tv-i, òn} < í=^ > {b®}- 



5.2 Ort honor mality of the {E® b }. 



The operators E® b form, with respect to the Hilbert-Schmidt norm, a set of d 2 orthogonal 
operators. The orthogonality of the j-E<® b j i s inherited from the orthogonality of the oper- 
ators Í-E'a^òfe}' which is shown in appendix A, equation (XIV). Using properties of tensors 
from M, we have 



^0,6' E g,h 



Trace 
Trace 
= Trace 
= Trace 



Kk ® e 



(2) 
0.2, í>2 



Trace 

(TV) \t 
ajv,fejv 



^o,6 ^s.h 



(X) 



01,61 02,02 



£(1) 

ai,bi 31,61 



(TV) 
(TV) 



02,62 92,"2 ' 



(7V)t 



i]V,6jv ÇN,hN 



E 



p(i) f p (i) 



Trace 



p(2)t F (2) 
02,62 gi,hi 



■ ■ ■ Trace 



ojv,6jv gN,h-N 



(dl ^01,31^61,^1) (^2 <5a 2 , 32^62, fc a ) "' (d'N 5 aN ,g N 0~b N ,h N ) — dS^gS^h, 

where we used our map between the sets {o^, 6 (fc) } -» {a®, 6®}, and the fact d = n*=f 4- 

íuAfc. The orthogonality of the 



^a,b>^g,h 



Thus we conclude (E® b , Eg h ) = Trace 
{ E a,b} means we can expand in terms of the {^® fe }, yielding p® = | Eo,6e{0,i,2,-,d-i} «0,6 #® 



Another property of the we shall need is the result of E® h E® b E® h , Using equation 

( jynjD , and the tensor nature of ^S» we have E f,h E a,b E fl = 



(2) 
02,62 



(8) • • • (g> -E. 



(TV) 
ojv,6jv 



e (1) , ®í; (2) , (8 

31, hi 92,6.2 ^ 



0£ 



(TV) 
9jv,6jv 
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ÍE {1 \ ®E {2 \ ®---®E {N \ )(eP\ ®E {2 \ ®---®E {N \ )(e {1) ! <8>E?>1 8- 

V 9lM 92,h 2 9N,h N ) \ ai,bí a 2 ,b 2 a N ,b N ) \ giM g 2 ,h 2 

E (1 \ E {1 \ É®1 ) ® ( E^ 2 \ E {2 \ eF>1 W • • ® ( E {N \ E {N \ E^l ' 

giM 01,61 giM I \ 32, fe a.2,02 32, fe I \ gN,n-N «jv.ojv 3jv,"jv 

= ^ c ® <!o 2 ® • • • ® <L = ^ c <6> (XI) 

where w fc = e d fc , = e 2 ?, and c = Efc=f ("fc^-fc - b k g k )^. 



5.3 The channel £ ® is unital and diagonal in the basis. 

The channel £® is diagonal in the J5® 6 basis. To see this, note that 

to) = *° k; 6i ® <k ® • • • ® <d (xii) 



= ® ® • • • ® £ W ^ ... ^ 

= («S,i <k) ® («2,2 * ■ ■ ■ ° KJL <L 

~ a axM a a 2 M " " " a «w,b W ^ai,6ij ® ^ a2 ,6 2 J ® • • • ® ^ ajV)6jv 

= a„, 6 «; 6l ) ® (*s 2 ) ® - ® «L) = A a , 6 ^ 6) 

where A Qj ft = a^ 6l 6 2 '"" a i^òjv' anc ^ we use d our bijective map ja^,ò^H <í=í> 
{a®,6®} to move back and forth between the operator basis set {-Efj,} and the opera- 

tor basis set {-E^* 6l ® E^ b2 ® • • • ® 4^ 6 N }' ^hus * ne tensor product of diagonal qudit 
channels yields a diagonal qudit channel. 

Next note that Ef = É$ ® ® • • • ® = I dl ® J d2 ® • • ■ ® I ajv = h- Taking a 



special case of the result in equation (|XIJ), we obtain 

£° (I d ) = £® (^o) = (42) ® ^ (2) (*§) ® • • • ® £ {N) «?) 

= (J dl ) ® £^ (I d2 ) ® • • • ® (/ djv ) = J rfl ® J d2 ® • • • ® I djv = / d . 

We conclude that £® (Id) = Id, and the channel £® is unital. Thus the tensor product of 
diagonal, unital qudit channels yields a diagonal unital qudit channel. 
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As an example, consider the product of two qubit (diagonal) unital channels, 8^ with 
diagonal parameters {Ai , A2 , A3}, and £^ with diagonal parameters {£1 , £2 , £3}- The 
product channel 8® = <g> £ (2) is a diagonal, unital channel, taking an input vector of 
(di (I2) 2 — 1 = 4 2 — 1 = 15 input density matrix coeficients a a>b to the output density 
matrix coefficients à a: b, as shown below. 



{ basis element I2 


X> a x } 


"«0,1" 




"«0,1 


= Cl «0,1 " 


{ basis element I2 




«0,2 




«0,2 


= 6 «0,2 


{ basis element I2 


®<T Z } 


«0,3 




«0,3 


= 6 «0,3 


{ basis element a x 




«1,0 




«1,0 


= Ai «1,0 


{ basis element a y 


®h] 


«2,0 




«2,0 


= A2 «2,0 


{ basis element a z 


®h } 


«3,0 




«3,0 


= A 3 a 3i0 


{ basis element a x 


<8> cr x } 


«1,1 




«1,1 


= Ai d 0:1,1 


{ basis element a x 




«1,2 


S ^ 


«1,2 


= AiÇ 2 «i,2 


{ basis element a x 


®a z } 


«1,3 




«1,3 


= AiÇ 3 ai i3 


{ basis element a y 


®a x } 


«2,1 




«2,1 


= A 2 Çi«2,i 


{ basis element a y 


®<J y } 


«2,2 




«2,2 


= A2 Í2 «2,2 


{ basis element a y 


<8> a z } 


«2,3 




«2,3 


= A2 Ç 3 «2,3 


{ basis element a z 


®a x } 


«3,1 




«3,1 


= A 3 Çia 3i i 


{ basis element a z 




«3,2 




«3,2 


= A 3 Ç2 «3,2 


{ basis element a z 


®a z } 


-«3,3- 




-«3,3 


= A 3 Ç 3 a 3j3 . 



5.4 The average output state of an optimal ensemble $ is oc Id for E®. 

Define the set of d 3 operators {F® bc } as F® bc = e^ E® b . Using our bijective map 
between {a±, 02, ci 3 , ■ ■ ■ , ajv-ij a N} and {a}, we expand F® b c in terms of a phase e~ and 
the {-Ea^}' ^ ur ex P ress i° n f° r F® bc becomes 

F® = e^E® = e™E? b ®E n 2 \ ® • • • ® E^\ . 



2iric 



The set of operators {F® b c } are the product of a phase e <i and the tensor products 

of the individual operators {-^a^b*,}- ^he {F® bc } are unitary operators, inheriting this 
behavior from the unitary nature of the phase factor and the unitary nature of the subsystem 
operators {-E^^j. To see this, note 

Ff bc F^ bc = (e^E a X \ ®E[ 2 \ ®-..®E a N \ V (e^ e£\ ® eP\ ® ■ ■ ■ ® E^\ ) 

a,b,c a,b,c y a\,bi 02,02 iJV,ojv / \ ai,í>i a2,»2 ajv,ojv / 

= =™ í E iiV 9E pV ®...®EÍ N V )(e£\ ®e[ 2 \ ®---®EÍ N \ ) 

I 01,61 a 2 ,b 2 a N ,b N J \ ai,í>i a 2 ,b 2 a N ,b N J 

= 1 í 1 ® ( E i2) l E M h )® ■ ■ ■ ® ( E W ! E^ N \ ' 

1 ai,6i ai,6i / \ a 2 ,b 2 a 2 ,b 2 I ^ ^ \ a N ,b N a N ,b N 
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= I dl <g> I d2 <g> ■ ■ • <g> I dN = Id, 

Í(k) ~| (fc)t (k) 

a k b f to say E a b E a b = . 

The {F® b c } form an irreducible group which we shall call Q. To see why Q is irreducible, 
recali our relation for irreducibility from j| discussed above. A necessary and sufncient 

1 2 

condition for a finite group <7 to be irreducible is if the relation jrL· J2 g eg Trace[g] = 1 
is trueQ. Here \\Q\\ is the order of the group Q. Let the group Q be the set {F® bc }, wnere 
a,b,c E {0, 1, 2, . . . , d — 1}. Q is of order d 3 and hence finite. Previously, we noted that 



Id and Trace 
the Trace sum yields 



o.o 



I7.®t rp® 



Saqftbh- Thus Trace 



E 



a,b 



dSaflSbfi. Computing 



oir E | T ™ce[ç] 2 = 5] 



\Q 



qeQ 



E 



E 



c6{0,l,2,...,d-l} 66{0,l,2,...,d-l} ae{0,l,2,...,d-l} 



Trace[F^ 



b.d 



- y 

d 3 z - 



E 



E 



ce{0,l,2,...,d-l} 6e{0,l,2,...,ci-l} ae{0,l,2,...,d-l} 



Trace [e E 



a.òJ 



^3 E E E 

c6{0,l,2,...,<i-l} 6e{0,l,2,...,d-l} oe{0,l,2,...,d-l} 
1 

" d3 



2iri 

e d Trace 



a,bi 



E E E 

cG{0,l,2,...,d-l} 6e{0,l,2,...,íi-l} ae{0,l,2,...,d-l} 
1 



Trace{E® ( 



d 3 



d 



E 



E 



6e{0,l,2,...,d-l} ae{0,l,2,...,d-l} 

Thus we find the group Q is irreducible. 



dd afi 5, 



b,0 



1. 



The fact that the channel £® is diagonal in the operator basis coupled with the 



equation (Q) result that E® h E® b E® h = ü c E® b , and the equation ( gg ) result that 
£ \E® b j = K a ^E® b , allows us to conclude the operators Í-P^bc} ano - * ne channel 



commute. 



*?M Z (P) F fL = E f,h £ (P) Ef h = Ef A ( i £ «a,6 A a , b j íJÍ (XIII) 

= 3 E a <^ Aa > b E 9,h E tb E gl = £ ( E g,h P E g,h) = £ ( F g,h,j P F fÍ,. 



Note that the product channel analysis in equation ( [XIII ) is essentially the same derivation 
as was done in equation (|ÏX| ) for qudits in the X a Z b operator basis. 



This is the key criterion for ensemble achievability. Since the ^F® b c | are unitary, F® h ■ pF® ] h ■ 
is a vàlid density operator. Applying any member of {F® b c } to an output optimal ensemble 
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{Pi > pf} yields an achievable ensemble. Since the group {Ff b c } is irreducible, we can apply 
Schur's lemma and conclude the average output state $® for an optimal ensemble for the 
product channel £® must equal -j Id- 

The remainder of our analysis for diagonal unital qudit channels uses the Schumacher and 
Westmoreland results summarized in equations (|ÏV|), (|y|), flVip and ( [VI 1| ) in the manner 
seen previously, and directly carries over to the product channel case. Thus we conclude 
for the product channel £®, the HSW channel capacity is 

N 

C = log 2 (d) - rmn p S(£ 9 (p)) = £ log 2 (d fc ) - min p S(£®(p)) . 

k=l 

6 Discussion and Conclusions 

The HSW channel capacity for single qubit unital channels was originally derived in ||] as 

C = 1 — miripS (£(p)) . 

This result was extended in |J to the tensor product of single qubit unital channels. For 
qubits, it was shown in || that there always exists a special basis in which a qubit unital 
channel can be written in diagonal form. A key step in their proof was a homomorphism 
between SU(d) and SO(d 2 — 1). Such a homomorphism would be necessary for the method 
of proof in H to carry through to the general qudit case for d > 2. However this homomor- 
phism only occurs for d = 2. Our method for deriving the HSW channel capacity depends 
on the qudit unital channel being diagonal, so our method only allows us to conclude that 

C = log 2 (d) - mín p S(£(p)) 

holds for diagonal unital channels. However, our proof was handcrafted in two key respects. 
The first was the choice of a fixed operator basis, the Generalized Pauli basis, in which the 
density matrix expansions were made. There exists the possibility that, given a specific 
channel, a custom operator basis could be constructed in which the channel £ would be 
diagonal. This in essence is how the proof showing any unital qubit channel is diagonal in 
some operator basis, was done in Q. 

The second assumption was the explicit manner by which we showed ensemble achievability. 
To summarize, we showed an output ensemble was achievable by 

1) restricting our attention to a preordained unitary operator basis consisting of elements 
and 

2) considering only diagonal channels in the basis Q. 
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The result was an algorithm by which we were able to determine, given an optimal ensemble 
ce, if the output ensemble was achievable for g G Q. 

The possibility remains that, given a channel £ , we could use a technique other than that 
developed in this paper to assure ensemble achievability across all elements of a group Q 
acting on the channel outputs of an optimal input ensemble. Again, this is essentially what 
occurs in the unital qubit channel scenario analyzed in Q. 

As a result, we feel we have "overconstrained" the requirements for our proofs. We conjec- 
ture the relation 

C = log 2 (d) - min p S(£(p)) 

holds for all unital qudit channels, rather than just those unital channels which are diagonal 
in the Generalized Pauli basis. 

As our final remark, the diagonal unital qudit channel capacity result extends the connection 
between the minimum von Neumann entropy at the channel output and the HSW channel 
capacity, which had previously been established in the qubit case, to a non-empty set of 
channels in any dimension. This implies a more universal connection between the minimum 
von Neumann entropy at the channel output and the classical information capacity for that 
quantum channel than had previously been shown. 

Furthermore, recali that the Holevo quantity x utilizes von Neumann entropy to obtain a 
relation for the distinguishability of quantum states. Hence it is reassuring that von Neu- 
mann entropy appears explicitly in our qudit channel capacity result. This is an indicator 
of consistency that reaffirms the fundamental role von Neumann entropy appears to play 
in Quantum Information Science. 



7 Acknowledgements 

The author would like to thank Patrick Hayden, Eric Rains and David Bacon for fruitful 
discussions during the development of this paper. 



17 



A The Generalized Pauli Group 



The generalized Pauli operators X and Z are used in our qudit analysis. This section 
describes some of the properties of these operators. Their definitions are 

X\j) = \j + l(modd)} and Z\j) = SV\j). 

The quantity = e~ . Note that X = Z = I^. The commutation relation of X 
and Z follows directly, yielding ZX = QXZ. Using the fact that (j + = 1, taking 

the Hermitian conjugate of both sides yields {j\X^\j + 1) = 1, allowing us to conclude 

= |i — 1 (modd)}. This in turn implies X is unitary, since XX^ = X^X = 
Similarly Z^\j) = Q^i \j), from which it follows that Z is a unitary operator. 

In our application of Schur's Lemma, we use the operator set of E a ^ = X a Z b , where 
{a, b} = 0, 1, 2, ■ ■ ■ , d — 1. We shall also use the operators F a ^ c = ü c X a Z b , where 
{a, b,c} = 0, 1, 2, • • • , d— 1. The operators E a ^ and F a ^^ c are unitary, since the composition 
of unitary operators is unitary. Note that E\ b = Z~ b X~ a and F^ bc = ÇL~ c E\ b . 



We now show that any qudit density operator p can be expanded as 
P=\ E c^X*& = \ E 

a,6e{0,l,2,-,d-l} a,&e{0,l,2,-,d-l} 



where the a a ^ are complex quantities. We shall work in the Hilbert-Schmidt operator norm, 
which for qudit operators A and B is defined as (A, B) = Trace[A* B\. Define the rescaled 
operators Q a ^ = -jà = x ^ . The operators Q a ^ are a set of d 2 orthonormal operators 
in the Hilbert-Schmidt inner product, as shown below. 

{Qa,b,Q q ,r) = - d {E atb ,E q , r ) = - d Trace\E\ h E^\ = ^Trace^X^ X" Z r ] (XIV) 

1 1 d_1 

(By the cyclic nature of trace) = -Trace[X q ~ a Z r ~ b \ = (j\X q ~ a Z r ~ b \j) 

d d j=0 

= 1 g n (r- b)j {j \x«-"\j) = I]T Çl^{j\j + q-a{modd)) 

3=0 3=0 



= 2^ a ' q E ^ j = -^d5 a , q 6b,r = à a ,q 6b,r- 

3=0 

Here 5 a ^ is the Kronecker delta function. Recali any qudit density operator p can be written 



as 

d-l d-l 



P = E E M a )( b l 

a=0 6=0 
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where the /3 ai b are complex quantities. We shall show that \a){b\ may be written as \a)(b\ = 
J2r=o Ss=o Cr,sQr,s, where the Cr,s are complex quantities. Rescaling the ( rtS , we will 
conclude that p may be written as 



d-l d-l 



P = E E a a,bE a ,, 



a=0 6=0 



To begin, write Q r ^ s as 

Define Cr,s as || 

( a h = Trace 



d-l 



i=0 



QÍJa)(b\ 



Vd 



Trace 



d-l 

E ^ s |i)(i + r|a}(6| 

3=0 



(XV) 



= ^Do i/ie Trace in the basis — * — /= E E ^ ^ (^|j)0 + r l a )(^h 

v« i=0 j =0 

-y d—1 d—1 ^ d—1 ^ 

where <5 is the Kronecker delta function. 

Consider the operator L = \a)(b\, and the corresponding complex coefficients Ç rjS 
(Q r ^s,L) = Trace Q\. s \a){b\ . We would like to expand L as L = J2r,s (Qr,s > L)Q r ^ s 

J2 r s ír,sQr,s- Note that ||L|| = 
conclude that 



-1 d-l 



(L,L) = 1. Using the result of equation ( |XV|) , we can 



E E \U 2 = E E 



i 

V^ 



d 



E E l*, 



6+r | 



E 



6+r | 



1. 



Thus 5Zr |Çr,s| 2 = 1 = ||L|| 2 . This fact for arbitrary a and b in |a)(ò| allows us to 
conclude the Q r)S form a complete, orthonormal basis for the L's, and we can expand L in 
terms of the Q r ^ s V o, 60]. Thus the expansion \a)(b\ = J2 r s {Qr,s , L) Q r ^ s holds V a, b. 
This leads to an expansion for the qudit density operator p. 

d-l d-l d-l d-l d-l d-l 

p = E E Mo)(b\ = E E P*,b E E (Qr,s,(\a)(b\) )Qr,s (xvi) 

a=0 6=0 a=0 6=0 r=0 s=0 

d-l d-l d-l d-l d-l d-l I /d-l d-l \ \ 

= E E E E 0a,b(Qr,s,(W)(b\) )Qr,s = E E ( Qr,s , E E ) Qr,s 

r=0 s=0 a=0 6=0 r=0 s=0 \ V a=0 6=0 / / 
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d-1 d-1 d-1 d-1 d-1 d-1 v -i d-1 d-1 

= E E (<w)q,, = E E 7f Or,. = E E 7f % = EE «r,.^,- 

r=0 s=0 r=0 s=0 V " r=0 s=0 V« V« " r= o s =0 

where — ^ = {Q r ,s , p) or equivalently a TyS = (E TyS ,p). 

The linearity of the inner product in the second argument was used to move the sum over 
the indices a and b inside the inner product. 

To obtain the final form of the expansion for the qudit operator p we shall use, note that 
Eq,o = Id- Om* result above, (E a ^ , E qj1 .) = Trace[E^ ab E q , r ] = dò as S^ r , telis us that 
Trace(E ab ) = dS a fl^b,o- Thus of the d 2 possible E a ^, only Eq$ has nonzero Trace. The 
trace condition Trace(p) = 1 allows us to conclude ao,o = 1- Using this, let T denote the 
set of d 2 — 1 elements a,b G {0, 1, 2, • • • , d — 1} with the exception that a and b cannot both 
be zero. Then we may write the qudit density matrix p as p = \ (i^ + J2( a ,b)er a a,b Ea,b) 
with a a>b = {E a)b , p) = Trace[E\ b p]. 

In the expansion of p above, there are 2d 2 — 2 real, independent degrees of freedom in 
the set of coefficients a ayb . However, in the density operator p, there are only d 2 — 1 real, 
independent degrees of freedom. Hence there are constraint relations between the a ajb . 

These constraints arise from the Hermitian nature of p. Note that E^ b 

Z-bj(-a = n d-b)(d-a) x d-a Z d~b = n d-b)(d-a) E d _ a4 _ b . Consideration of pt = p then 

implies 





= \[h+ E < b ^ d - a)(d - b) Ed-a,d- b ) 

V (a,b)er J 

or ad- a ,d-b = ot* a b íí( d_a )( d_6 ) . Here * indicates complex conjugation, and index arithmetic 
is modulo d. 

For example, for qubits, d = 2, and Q = e~ = e m = —1. Applying the constraint 
equation above leads to a* 1 ^( 2 -°)( 2 - 1 ) 

= «2-0,2-1 or oíq i = ao,i) implying the coefficient 
of £?o,i = Z must be real. Similarly, o^gí^ 2 " 1 ^ 2-0 ) = «2-1,2-0 or a* = 0:1,0, implying 
the coefficient of E\$ = X must be real. Lastly, aj^í^ 2-1 ^ 2-1 ) = «2-1,2-1 ° r ~ a i y i = 
ai,i, implying the coefficient of E\ y i = XZ must be pure imaginary. Note that X = a x , 
XZ = —i<J y , and Z = a z . Hence we have reproduced the Bloch Sphere representation for 
qubits, p = \ (l 2 + a lfi X + cxi^XZ + cto.i-Z) = \ (h + w x a x + iw y (-i<r y ) + w z a z ), 
with the Wk real. For qubits, we end up with 3 = d 2 — 1 real independent parameters, 
and not 2cí 2 — 2 = 6. The constraint equations for the a ayb eliminated three real degrees 
of freedom. In general, the constraint equations will eliminate d 2 — 1 real extra degrees of 
freedom, leaving d 2 — 1 actual real parameters. 
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